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Abstract 

We study critical and universal behaviors of unitary invariant non-gaussian ran- 
dom matrix ensembles within the framework of the large-iV renormalization 
group. For a simple double-well model we find an unstable fixed point and a 
stable inverse-gaussian fixed point. The former is identified as the critical point 
of single/double-arc phase transition with a discontinuity of the third derivative 
of the free energy. The latter signifies a novel universality of large- N correlators 
other than the usual single arc type. This phase structure is consistent with 
the universality classification of two-level correlators for multiple-arc models by 
Ambj0rn and Akemann. We also establish the stability of the gaussian fixed 
point in the multi-coupling model. 
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1 Introduction 



The universal behaviors in random matrix theories are key ingredients in their 
level-statistical application, and have been discussed extensively. It is known 
that there are several universal quantities which do not explicitly depend on the 
probability distribution of a matrix ensemble. Because of the universality we may 
calculate such quantities from a simple gaussian ensemble which shares the same 
symmetries as the physical system in concern. Recently, Ambj0rn and Akemann 
pointed out a universality classification of matrix ensembles with respect to the 
two-level correlators []J. The large- N (or "smoothed" ||) connected two- level 
correlator depends solely on the support configuration of the level density p(A). 
For example, this correlator in a matrix ensemble with a single-arc level density 
is identical to that calculated in the gaussian ensemble. 

In this letter, we study the phase structure of a non-gaussian matrix ensemble 
with a multiple-arc level density in terms of the large- N renormalization group 
[Q. The general algorithm of the large- N renormalization group enables us to 
find out the exact location of the critical point as a fixed point and to specify 
the order of the phase transition exactly [|J . Moreover, its renormalization group 
flow clarifies the global phase structure by classifying the parameter space into 
phases separated by critical surfaces emanating from unstable fixed points ||. 
Along this line, a conjectural account for the branched-polymeric behavior of 
c > 1 bosonic strings was recently provided || . 

We first concern ourselves on the merging phase transition where two energy 
bands coalesce into a single band. This merging transition has recently attracted 
attention, for instance, in the application of random matrix theory to the level 
statistics of QCD Dirac operator 0; there the level density around A = is 
identified with the chiral condensate of quarks, and thus its vanishing is related 
to the chiral symmetry restoration at a finite temperature |J. More specifi- 
cally, in order to study the phase structure involving different arc configurations, 
we investigate the simplest non-gaussian model with a U(iV) and Z 2 -invariant 
distribution P(4>) oc exp — iVtr (— 10 2 + |0 4 ) , which can take either single or 
double-arc level density depending on the parameter g. It is known that at a 
merging point g = 1/4, the third derivative of the free energy has a discontinuity 
due to the switchover of the functional form of p(A) ||. We shall correctly identify 
this critical point as an unstable fixed point of the renormalization group, from 
which the renormalization group flow is emitted. As a straightforward exercise 
we derive the explicit form of the free energy by integrating the renormalization 
group equation. 

Next we shall discuss the off-critical, universal behavior of the model. It is 
argued to be determined by the nature around the stable fixed point which is a 
final destination of the renormalization group flow 0. To proceed we need to 
define the approximate /3-function with required properties. The phase diagram 
described by such a /3-function is in accord with the above mentioned universality 



1 



classification by the large- N two-level correlator. Finally we establish the stability 
of the gaussian point by calculating all the associated scaling exponents in multi- 
coupling models. 



2 Renormalization group approach 

We consider an ensemble of random N x N hermitian matrices <fi with a weight 



P N {4>) = Z^rW-^trV^)], V{<j>) = £ f <p k . 

k=2 K 



(1) 



Hereafter the coupling constant g<i is fixed either to —1 (in sect. 3, 4, 5) or to +1 
(in sect. 6), and (72-dependences will not be indicated explicitly. The expectation 
value with respect to the weight is denoted by (••■). The partition function 
Zpf(gj) is defined so that (1) = 1. It is written as an integral over eigenvalues 



f N 2 
z N{gj) = c N Y[d\ e Y[ - Aj| 2 exp 



l<i<j<N 



N 



k=l 



(2) 



after an integration over angular part of <p, which leaves us the volume cat of 
U(iV). The free energy is defined as 



ZN{gj 



ZN(gj 



IS2 = lJ 



(3) 



We briefly recall the renormalization group approach to random matrix mod- 
els developed in refs.p, £|, |5|. In the spirit of the renormalization group, we 
integrate a part of the degrees of freedom, namely the (N + l)-th eigenvalue 
Ajv+i = A, in Z^+iigj) in order to relate it to Z^(gj + Sgj), 



ZN+\{gj 



CjV+l 

cn 



(N+l)trV{<f>) / ^ e -(7V+l)y(A)+2trlog|A- 



CN 



Z N ( 9j )exp [-(tr Vm - W«A S » + 2(tr log |(A S ) - <f>\)\ . (4) 



In the second line, we have evaluated the A-integration by the large-A^ saddle 
point, and used the factorization property of the correlation functions. It leads 
to the saddle point equation for the expectation value (A s (#.,•, </>)): 



n<A s » 



tr 



N (A s )-0 ( 



(5) 



By the use of the loop (Schwinger-Dyson) equation, eq.fl5|) is shown to be equiv- 
alent to the condition 



P«As» = 0, 



p(A) := — — Im ( — tr 



1 



N X + ie-(p l 



(6) 
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which means that the (JV+l)-th eigenvalue settles down to one of the edges of the 
arcs of p(A). Expanding eq.(|J) into a series in 1/N, we have the renormalization 
group equation obeyed by the free energy: 



d 

N 8N +2 



\ + trV(^)) + V((X S )) - 2 tr log |(A S ) - 0|) + O (i 



G ft -]- (7) 



<9F* 
<%, 

In the right hand side the coupling constants g 3 , . . . ,g m and the one-point func- 
tions a,j '■= §g~ = j( ^ tr J ) (3 < j < m) are regarded as independent variables. 
Other one-point functions appearing in the right hand side are understood as 
functions of g 3 - and aj (3 < j < m) via the Schwinger- Dyson equations. In the 
same way, (A s ) is also regarded as a function of gj, aj (3 < j < m). The explicit 
form of G is given by 

G( 9j , a s ) = -1 + f; (l - J -) g,a 3 + V((X S )) (8) 

- 2 C d2 ((h i 7 L i)-\)- 2, °^ 

Let us turn to the process of extracting singular behavior of F from the 
renormalization group equation (|7|). For simplicity, we consider the case with 
a single coupling constant g. The extension to multi-coupling constant case is 
straightforward ||]. We assume that the free energy in the large- N limit has a 
fractional power singularity 

oo oo 

F(g,N) = Y,a k (g-g*) k + J2h(g-g*) k+J + 0(N- 2 ) ( 7 e R+\Z, b ^ 0). (9) 

fc=0 fc=0 

The critical point g*, critical exponent 7 and coefficients a±, a 2 are determined by 
solving the following closed set of equations 

= G A (g m ,ai) , (10a) 
2 

- = G ga {g^ai) + 2a 2 G, a>a (g*,a 1 ) , (10b) 
7 

2^ = G, g {g*,ai) , (10c) 

2a 2 = -G i9i9 (^,ai) + 2a 2 G i9ia (5f*,ai) + 2(a2) 2 G iaia (^,ai) , (lOd) 

where j9 and >a denote partial derivatives. These condition can easily be estab- 
lished by substituting (|) into the large- iV renormalization group equation 

2F(N,g)=G(g,^Y (11) 
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By comparing coefficients of {g — g*) k and (g — g*) h+1 in the both sides we obtain 
a family of equations for aj, bj, g*, and 7. The coefficients of (g — g*) k+1 contain 
both a,j and bj, whereas those of (g — g*) k do only a,j (thus they can be solved 
for a,j). Since these equations are homogeneous in bo the overall normalization of 
the singular part is left undetermined. 

The above procedure can as well be used to detect the discontinuity of a 
derivative of the free energy. We now assume that the 7-th derivative of the free 
energy is discontinuous: 

00 00 

F( g , N) = J2 a *(9 - g*) k + E - g*) k+1 ^[g - 9*} + 0(N- 2 ) (12) 

k=0 k=0 

( 7 eZ + ,6 ^0). 

By substituting (|T2"D into the large- N renormalization group equation and com- 
paring coefficients of (g — g*) k and (g — g*) k sgn[g — g*], we can easily show that 
eqs.(10) holdE. 



3 Merging transition fixed point 

The simplest model which undergoes the merging transition is the one with quar- 
tic potential with negative quadratic term |§ 

V{<P) = -\<P 2 + \<P'. (13) 

Though a parameter g 3 is missing compared to (0), we can set g Q dd = a dd = 
consistently because of the Z2 symmetry^. 

We present the explicit form of the G-function: 

G(g, a ) = -l-ga- i(A s ) 2 + |(A S ) 4 - 2 log(A s ) (14) 

r(>*) ( 1 2\ 

-J dz i-z + gz 3 - yj(-z + gz 3 ) 2 - 4(-l + gz 2 - g + ^g 2 a) J 

which is supplemented by the saddle point equation: 

(-(A s ) + s(A s > 3 ) 2 - 4(-l + ^(A s ) 2 - g + 4^ 2 a) = (15) 
determining (A s ) function of g and a. 



*So far we have implicitly assumed that the powers of the singular parts do not contribute 
to the regular parts. It is true for 7 ^ Q or 7 = p/q with p — q > 3 and p, q are co-prime 
(fractional power singularity) and for 7 > 3 (discontinuity) . 

tHere we have excluded the possibility of spontaneous breaking of the Z2 symmetry 11 . The 
linearized RG flow for the broken Z2 phase is readily found in Fig. 3 of ref.|| by restricting it 
onto the curve 2g| — 9.94 = 0. 
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One can show analytically that g* = 1/4,7 = 3, a± = 5, 02 = —18 solves 
eqs.(10). The integer value of the 7 exponent suggests that around g* = 1/4 
the free energy is of the latter form ([12]) and has a discontinuity in its third 
derivative.. This is in accord with the large- N solution found in |§; at g \ 1/4, 
the double-well potential becomes well separated so that the level density splits 
into two disconnected arcs. We stress that our renormalization group equation 
can actually detect not only the fractional power behavior of F(g) but also its 
discontinuity due to the switchover of the functional form of p(A). 

The readers should notice that, in the previous procedure we have not made 
use of any kinds of multiple-arc ansatz, which is used in obtaining large- N solu- 



tions [|I(J. If we were interested only in a particular arc configuration, we could 
assume a factorized form of the discriminant of the square root in the G-function 
(|T4D. However our method, which uses the discriminant containing a variable 
a as it is, can represent any possible arc configurations. Thus by exploring all 
solutions to eqs.(10) in the two-parameter space (g, a), we can find solutions with 
arbitrary arc configurations respecting the Z 2 symmetry. 



4 Free energy 

To obtain the free energy, one usually reads off the level density p(A) from the 
expectation value of the resolvent and then takes the average of V(\) — V log |A — 
A' I with respect to it ||10|| . Here we show that the explicit form of the free energy 
can as well be obtained by integrating the renormalization group equation (0). 
Our calculation avoids computing cumbersome principal-part integration. 

For this practical use, it suffices to choose a particular arc configuration from 
the outset and to use the fact that the saddle point (A s ) falls on an edge of the 
arcs due to the level repulsion. Then (A s ) is expressed as a function of g only and 
the renormalization group equation takes a simple linear form. 

In g > g* case, assuming a single arc configuration at [—A, A] and the asymp- 
totics iV _1 (tr (z — </>) _1 ) ~ 1/z (\z\ — > 00), we obtain 




-z + gz 3 - (-1 + 9 -A 2 + gz 2 ^j Vz 2 - A 2 



(16) 



where A 2 = ^(1 + ^/l + 12g). Knowing that the saddle point is given by the 
edge of the arc (A s ) = ±A, the G-function (0) then takes the form: 

/ dF\ dF 1 A 2 , A 2 



In the leading order of 1/N, the renormalization group equation ([?]) becomes 

dF 
dg 



2F(g) = G(g,±P). (U 
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The solution to this differential equation with the G-function (O) is 

3 5A 2 A 4 1 A 2 d 
F ( 9 » = -8-^- 384 - 2 log T + F' < 19 > 

where C\ is an integration constant. 

In g < g* case, we assume the double arc configuration at [— A 2 , —A{\ U 
[Ai,A 2 ], 



1 1 \ 1 
tr 



(-z + gz 3 - gzy/(z*-A$)(z*-J%?) , (20) 



\N z-<j>l 2 

where A\ = (1 — 2^fg)j g and A 2 2 = (1 + 2^/g)/g. By requiring the continuity 
from the g > g* case, we are obliged to choose the outer edge ±A 2 as the relevant 
saddle point (A s ). Then the G-function ( 14] ) takes the form: 

_/ dF\ dF 1 1 1 



The solution to the renormalization group equation (|l~8|) with this G-function 

becomes 

F( 9 ) = -5 + il„ gff -i + ^ (22) 

with an integration constant C 2 . The integration constants C\ and G2 must be 
the same because of the continuity of F(g) at g = g*. On the other hand, C 2 must 
vanish since the free energy should be dominated by the minima of the potential 
in the small-g limit, 

F(s) ~ v ) = ~h (9X0) - (23) 

This result ( |T5[ ) and ( P2| ) with Gi = G2 = is identical to that obtained in ref . ||i 
using BIPZ method. 



5 Inverse-gaussian fixed point 

The concept of flow is indispensable in the renormalization group analysis. How- 
ever, the nonlinearity of our exact renormalization group equation (0) with re- 
spect to dF/dg makes an idea of flow ambiguous. Thus we wish to define an 
approximate flow which visualizes the phase structure correctly. Indeed there is 
such a sensible definition of flow. 

We define an approximate /3-function expanded around g as 

9G(g, a) 



da 



(24) 

o=ai+2a 2 (g-g) 



"'"Its free energies err: Eb misses a constant \ log 2; Eq should read (in our convention (|l3|)) 



3 , -l + ya-153 , -33-10 yi-15 3 1 i--/ l-yi-15 g \ 
'8 + 675^ 1 180^ 2 l °S{ —g )■ 
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where a± 

Let us first note that 



«2 



1 £F 
2 



(g), and shall check the validity of this definition. 

8G(g,a) 



(3 e ™ ct (g) 



da 



(25) 



is the 'exact' 
equation 



beta function in the sense that the free energy F obeys a linear 



N 



_d_ 

dN 



F(N,g) = r(g). 



(26) 



The above 'exact' beta function is practically useless because it requires the 
full information of the free energy. It is desired to extract correct qualitative 
information (the phase structure of the model or the universality classification 
of phase transitions) by the use of the renormalization group equation (|7]) in 
case the exact form of the free energy is not available. On the other hand, the 



approximated /5-function ( |24"|) is obtained by replacing dF/ dg with its truncated 
expansion around g. We can determine the coefficients a\ and ci2 by solving (10) 
even if full information of F is not available. Thus we can calculate /3(g) for 
arbitrary g only if eqs.(10) are solved at a single point g. 

Since we have kept the linear term in g in the power series expansion of 
dF/dg, the zeroth and the first derivatives of G at g = g are exactly reproduced. 
Therefore [3(g) enjoys the following propertyi: If g is a fixed point of the flow 
(i.e. solves eqs. (10)), then 



m 

d/3 





2 
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(27a) 
(27b) 



Even if g does not solve the equations (10), the definition (pl[) makes sense and 
f3(g) approximates the character of the renormalization group transformation in 
the vicinity of g. However, expansion around such a g is not very useful because 
of the lack of algorithm to determine a\(g). 

The definition (f?§) is easily extended to multi-coupling 



Pk(gj 



dG(g 



3 1 a jj 



da h 



- 2 Y.i a i^9l-gi) 



(28) 



where dj = dj(gk)- Notations a,j and aje are those used in ref.||. 

In Fig.|l], the /3-function for the model ( |13|) expanded around the merging fixed 
point g* = 1/4 is plotted. We observe that the approximate /5-function is smooth 



§ Previously we have proposed in Q a linearized /3-function 

/3 lincar (ff) := in 
oa 

Though this beta function vanishes at a fixed point g, its first derivative does not give the exact 



value of the critical exponent because of the missing second term in (10b) 
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Figure 1: Approximate /3-function (expanded around = 1/4). 



everywhere even though the matrix integral is ill-defined for g < 0. It vanishes 
at the very vicinity of the origin with a negative exponent. This indicates that 
the origin is an attractive fixed point. The level density becomes two infinitely 
separated semi-circles when approaching this point. 

It is argued by Brezin and Zee [§], although not rigorously proven, that the 
universality of the large- iV connected two-level correlator be attributed to the 
presence of attractive gaussian fixed point of the renormalization group transfor- 
mation. If we accept this argument, then the attractive inverse-gaussian fixed 
point in the model (|T3|) plays the same role for the new universality for the 
multiple-arc matrix model recently found in ref. 




6 Stability of the gaussian fixed point 

To strengthen our assertion that each fixed point corresponds to a universality 
class of large- N two-level correlators, we present here the analysis of stability of 
the gaussian potential in the mo del 

oo 

^) = Et^ 92 = +l. (29) 

k=2 K 

renormalization group analysis of the stability of the gaussian potential has been done 
by Y. Morita et. al. (Phys. Rev. B52 (1995) 4716) based on mapping from matrix models to 
vector models. They have imposed an arbitrary restriction to a smaller subspace of coupling 
constants which is inconsistent with the renormalization group flow, yielding incorrect values 
of the exponents. 
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We evaluate the derivatives of the G-function (§) at the gaussian point g k = 
(k — 3, 4, • • •). We immediately find that 



_d_ 

da t 



(As) 




9i=0 



implying 



0, a„) = G t a e ,a k (gj = 0, a n ) = 



(30) 



for arbitrary a n . Thus the gaussian point is always a fixed point. 

To obtain the exponents on the gaussian fixed point, we need to calculate 

contributes (1 — k/2)5e k to it. It can also be 

1 



G )9fciClr The second term in 
shown that 



d 2 



dg k da e 




Therefore the scaling exponent matrix G^ k : 
gular with eigenvalues 



if £ + 2>k. 



G ,9k,a e + 2a krG iae:ar is upper trian- 



G kk {0,0,...) 



3,4, 



(31) 



The negativity of the exponents indicates the stability of the gaussian fixed point 
in the entire coupling constant space (^) in accord with the universality of the 
large- N two-level correlator. 

In the (72 = — 1 case, the model at the inverse-gaussian point is ill-defined 
and therefore the general algorithm cannot extract any critical exponents around 
this point. Nevertheless, the model is well-defined except at the inverse-gaussian 
point. The inverse-gaussian point is a fixed point in the sense of the limiting 
procedure 

limG a (<?,a) = 0. (32) 

The stability of this fixed point is guaranteed by the negative definiteness of the 
extended /5-function G a for < g < |. 



7 Summary 

This letter is devoted mainly to the application of the method of large- iV renor- 
malization group to a matrix model with a double-well potential. The advantage 
of our method is that the arc configuration need not be specified from the outset, 
but rather a single renormalization group equation (|i~4T ) (and its approximated 
/3-function) encompasses all possible arc configurations. As a direct consequence 
of it, we have extracted exact values of the critical (fixed) point and exponent 
associated with the splitting/merging of the arc(s). Furthermore, by the use of 
approximate /5-function we have found a novel inverse-gaussian fixed point, whose 
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stability may guarantee the large- N universality of double- arc matrix ensemble. 
We have also calculated the scaling exponents around the gaussian fixed point 
for a matrix model with a generic potential. Their negativity, i.e. the stability 
of the gaussian point against perturbation of potential preserving the single-arc 
configuration, supports the large- N universality of (usual single arc) matrix en- 
sembles. 

We wish to utilize our results to the random matrix analysis of disordered 
physical systems. 
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